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Abstract 

The systematical studies on the dynamical approach of wavefunction collapse in quantum mea- 
surement are reported in this paper based on the Hepp-Coleman's model and its generalizations. 
Under certain physically reasonable conditions, which are easily satisfied by the practical problems, 
it is shown that the off-diagonal elements of the reduced density matrix vanish in quantum me- 
chanical evolution process in the macroscopic limit with a very large particle number N. Various 
examples with detector made up of oscillators of different spectrum distribution are used to illustrate 
this observations . With the two-level system as an explicit illustration, the quantum information 
entropy is exactly obtained to quantitatively describe the degree of decoherence for the so-called 
partial coherence caused by detector. The entropy for the case with many levels is computed based 
on perturbation method in the limits with very large and very small N. As an application of this 
general approach for quantum measurement, a dynamical realization of the quantum Zeno effect 
are present to analyse its recent testing experiment in connection with a description of transition 
in quantum information entropy. Finally, the Cini's model for the correlation between the states 
of the measured system and the detector is generalized for the case with many energy-level. It is 
shown that this generalization can also be invoked to give the dynamical realization of wavefunction 
collapse. 
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1. Introduction 

Though quantum mechanics has been experimentally proved to be a quite successful theory, 
its interpretation is still an important problem that physicists should face [1-4]. To interpret the 
physical meaning of its mathematical formalism, one has to invoke the wave packet collapse(WPC) 
(or wavefunction collapse)postulate as an extra assumption added to the closed system of laws in 
quantum mechanics . This postulate is also called von Neumann's projection rule or wavefunction 
reduction process. Let us now describe it briefly. It is well known in quantum physics that, if 
measured quantum system S is in a state \<f> > that is a linear superposition of the eigenstates \k > 
of the operator A of an observable A just before a measurement, ie., 

\<fi >— Ck\k >, c' k s are complex numbers (1.1) 

then a result of th measurement of A is one a k of the eigenvalues of A corresponding to \k > with 
the probability |cfc| 2 . The von Neumann's postulate tell us that, once a well-determined result 
about A has been obtained, the state of S is no longer \<f> > and it must collapses into \n > since 
the immediately-successive measurement of A after the first one should repeats the same result. In 
terms of the density matrix 

p=\ ( p>< ( t ) \=J2ckC* k ,\k><k'\, (1.2) 

k,k' 

for the state \<f>>, the above WPC process can be expressed as a projection or reduction 

/ o^P = ^|c„| 2 |fc><fc|. (1.3) 

Because the off-diagonal elements represent coherence, through which the density matrix describes 
a non-classical probability, the wavefunction collapse characterized by vanishing of off-diagonal 
elements means the loss of quantum coherence or called quantum decoherence [5] . 

The recent studies on the quantum decoherence in an open system S surrounded by an envi- 
ronment E was mainly motivated by the interests in the macroscopically-quantum effects such as 
dissipation in the quantum tunnelling and the semiclassical gravity theory for particle creation in 
cosmology [6-10]. In quantum dissipation theory, an important treatment for quantum decoher- 
ence is invoking quantum Brownian motion through its master equation to describe the vanishing 
of the off-diagonal elements of reduced density matrix of 5* [1]. The recent investigations in this 
context were carried out with the Feynman-Vernons integral method for the Ohmic, sub-Ohmic 
or super-Ohmic environment (e.g, see ref.[8]). Without the use of path-integral [9,10], an exactly- 
solvable dynamical model of quantum dissipation was presented by Yu and this author to deal with 
the similar phenomenon. Zurek and his collaborators especially emphasized the role of environment 
surrounding the open system, which monitors the observables of the system so that their eigenstates 
continuously decohere and then approach classical states [4,8]. 
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It should be noticed that there exists another dynamical theory based on Hepp-Colemen's inves- 
tigation to dynamically realize quantum decoherence characterized by wavefunction collapse [11,12]. 
This theory and its generalization [13-15] were proceeded in a purely-quantum mechanical frame- 
work. In fact, to realize the WPC, the external classical measuring apparatus detector must be 
used to detect the result. Then, someone thinks the WPC postulate to be not quite satisfactory 
since quantum mechanics is expected to be an universal theory valid for whole 'universe' because 
the detector, as a part of the universe, behaves classically in the von Neumann's postulate. A 
reasonable description of the detector should be quantum essentially and it exhibits the classical or 
macroscopic features in certain limits. If one deal with the detector as a subsystem of the closed 
system (the universe C = the measured system S + the detector D), it is possible that the quantum 
dynamics of the universe can result in the WPC through the interactions between S and D. Up 
to new, some exactly-solvable models have been presented to analyse this problem. Among them, 
the Happ-Coleman (HC) model is very famous one and has been extensively studied in last twenty 
years [13-19]. In order to describe studies in this paper clearly, we need to see some details of this 
model. 

In the original HC model, an ultrarelativistic particle is referred to the measured system S while 
a one-dimensional array of scatterers with spin- 1/2 to the detector D. The interaction between S 
and D is represented by an homogeneous coupling 

N 

H I = Y J V(x - a n )a[ n) (1.4) 

71=1 

(n) 

where a\ is the first component of Puli matrix; a n is the position of the scatterer assigned to the 
n'th site in the array. The Hamiltonian for D is 

H s = cP (1.5) 

where c, P and x are the light speed, the momentum and coordinate operators respectively for S. 
This model is quite simple, but it can be exactly solved to produce a deep insight on the dynamical 
description of the quantum measurement process. Starting with the initial state 

1-0(0) >=5^c fc |fc>®|£>> (1.6) 

where \D > is pure state of D (it is usually taken to be ground state), the evolution state \ip(t) > 
for the universe C=S+D is defined by the exact solution to this model. Then, the reduced density 
matrix 

Ps (t)=Tr D (\m><m\) (1-7) 
of the measured system is obtained by taking the trace of the density matrix 

P (t) = \m >< m\ (i.8) 
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of the universe to the variables of D. Obviously, p s {t) depends on the particle number N of D. When 
N — > oo, i.e., in the macroscopic limit , p s (t) — ► p after long enough time t as eq.(l-2). Namely, the 
Schrodinger evolution of the universe C=S+D leads to the WPC for the measured system. More 
recently, the original CH model was improved to describe the energy exchange between S and D by 
adding a free energy Hamiltonian 

N 

H = hujJ2 (T 3 n) ( L9 ) 

n=l 

and correspondingly improving the interaction slightly. Notice that the improved model remains 
exactly-solvable [13]. 

Because the spin quantum number is fixed to be 1/2 in the original CH model or its improved 
versions, they can not describe the classical characters of the measurement. Usually, the classi- 
cal feature of a quantum object is determined by taking certain value for some internal quantum 
numbers of the detector D or h = 0. In the case of the angular momentum, this classical limit 
corresponds to infinite spin. This generalized dynamical model was successfully built by this author 
in 1993 [14, 15]. The first step is to establish such a generalization of the HC model manifesting 
the WPC as the dynamical process in the classical limit as well as in the the macroscopic limit 
simultaneously Then, the essence for this model substantially resulting in the realization of the 
WPC as a quantum dynamical process as well as for those well-established was found to be the 
factorization of the evolution matrix in the interaction picture with help of a detailed study on the 
dynamics of the generalized HC model in both the exactly-solvable case and the non-solvable case, 
For the latter, the high-order adiabatic approximation (HOAA) method [20-23] is applied to its 
special case that the coupling parameter depends on the position of the measured ultrarelativistic 
particle quite slightly. Finally, we point out that this possible essence in the dynamical realization 
of the WPC, is largely independent of the concrete forms of model Hamiltonians. Notice that, in 
the dynamical models of wavefunction collapse for quantum measurement, both the macroscopic 
and classical measuring apparatus can be regarded as the environment in the certain model of the 
quantum decoherence for quantum dissipation [4,8]. This is because both they act as the classical 
or macroscopic monitors seeing the system. 

However, because all of the previous dynamical models of decoherence for quantum measurement 
depend on the specific forms of Hamiltonians of D and S, it is necessary to present a model- 
independent dynamical approach for decoherence in quantum measurement process based on the 
HC model. It is expected that such an approach does not depend on the detailed structures of the 
Hamiltonians of S and D as fully as possible, but can be invoked to deal with the practical problems 
of quantum measurement such as quantum Zeno effect (QZE). This universal approach should also 
be used to describe the role of environment in decoherence for quantum dissipation. It is shown in 
present study that, through a suitable choice of the interaction between S and D, the Schrodinger 
evolution of the universe C formed by S plus S may result in the phenomenon of decoherence in the 



4 



reduced density matrix p of S at the macroscopic limit of D with very large N. Mathematically, 
the mechanism of this phenomenon is that the accompanying factors in the off-diagonal elements 
of reduced density matrix p s caused by the dynamical evolution of C will vanish as N approaches 
infinity. It even was described in concrete examples [14,15]. Notice again that in the previous models 
for the quantum measurement and decoherence, the considered systems S usually are specified as 
an ultrarelativistic particle or a two-level system while the detector D as a spin array. Here, what 
we require is only that the system is of the non-degenerate discrete spectrum and the interaction 
between D and S is chosen to result in a factorizable evolution matrix. However, D is required 
to satisfy a condition that any row or column of each factor corresponding to the factorizable 
evolution matrix at least has one non-vanishing off-diagonal element. This means that the back 
action of the measured system can effectively act on the detector so that the microscopic states 
can be read out from the macroscopic counting numbers contributed by all particles in detector. 
This condition is physically reasonable and can be satisfied in widespread circumstances. Some 
examples are invoked to illustrate that this condition can be realized by choosing suitable spectrum 
distributions of oscillators making up the detector. 

To quantitatively describe the intermediate state of decoherence between the pure state and the 
most-largely mixed state, we need to calculate the entropy of S 

s = -yTr(pln / 3) (10) 

for two cases: (I) For the two-level system and any finite N, an exact solution for s is obtained as 
the functional of p. (II) For the limits with very large and very small N, the approximate solutions 
of s by certain perturbation methods. These calculations show that the entropies indeed decrease 
as N increases and they will take the maximum values at the infinite N. This means that an 
ideal macroscopic detector or environment must cause the increment of entropy of the its monitored 
system . 

The above general approach is used to built an exactly-solvable dynamical model for the quantum 
Zeno effect [25,26] in connection with the recent experiment by Itano e£.aZ[26] about the inhibition 
of quantum transition between the atomic energy levels. The present investigation compromises 
the different points of view about this experiment testing quantum Zeno effect [27-37]. In this 
model the detector is simplified as a system of N oscillators with a suitable interaction with the 
measured system-a two-level atom, we show that, due to gradually-vanishing of the off-diagonal 
elements in p$, the two- level system will be frozen in its initial level as the times of measurement 
in a given time interval becomes infinite. This is just the quantum dynamical realization of the 
QZE through a dynamical approach of the wavefunction collapse. The information entropy for the 
process of quantum Zeno effect of two-level system is calculated to manifest an interesting behavior 
of transition from random to regularity: for a given time interval, when the times L of measurement 
is less than a critical value L c , the entropy changes at quite random as L changes; when L is larger 
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than L c , the entropy decreases monotonically as L becomes larger. 

Finally, it has to be pointed out that the correlation between the states of the measured system 
and that of the detector has not been emphasized well in the original HC model and its generaliza- 
tion . But this problem was well analysed by Cini in his beautiful dynamical model [39] . The present 
investigation is also to emphasize on both the wavefunction collapse and the state correlation in 
the Cini model. In fact, the correlation between the states of measured system and the detector is 
crucial for a realistic process of measurement, which enjoys a scheme using the macroscopic counting 
number of the measuring instrument-detector to manifest the microscopic state of the measured sys- 
tem. Notice that the original Cini model for the correlation between the states of measured system 
S and the measuring instrument-detector D is build only for a two-level system interacting with the 
detector D, which consists of N indistinguishable particle with two possible states iv and u>i. For 
the two states u + and u_, the detector has different strengths of interaction with them. Then, the 
large number N of " ionized" particle in the ionized state uj\ transiting from the un- ionized state wo 
shows this correlations. In this paper, we wish to generalize the Cini's model for the M-level system. 

2. Dynamical Description of the General Model 

In this section we wish to describe a general dynamical model for quantum decoherence caused 
by a suitable interaction between the measured system S and a measuring instrument- detector (or 
an environment) D that can be regarded as a reservoir at temperature T. The considered system is 
only required to be of the non-degenerate discrete spectrum. Let \n)(n = 1, 2, ...M) be the discrete 
eigenstates of S corresponding to N energy levels E n (n = 1, 2, ...M). Therefore, the Hamiltonian 
is formally expressed as 

M 

H s = ^E n \n)(n\. 

n=l 

D is made up of ./V particles with the single particle Hamiltonian Hk(xk) for dynamical variables 
Xk (such as canonical coordinate , momentum and spin ) of the k'th particle. Its Hamiltonian 

N 

H D = Y J Hk(^k) = X>a|#fe(a*MaMa><0a| (2-2) 
k k.a 

can be written in terms of Hk(%k) or its eigenstates. Here, it has been assumed that there are not 
mutual interactions among detector's particles. 

Physically, the interaction between S and D can be chosen to have the different strengths for 
the different states of S. Thus, one can write the interacting Hamiltonian as 

H I = J2J29( t ) V n(^)\n)(n\ (2.3) 

n k 

where g(t) (= 1 for 0<i<r; = for t < or t > t) is a switching function; V n ^ V m for m =/= n. 
Besides the above-mentioned basic requirements, their is only a few of constrains on the model 
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and one even need not to know the specific forms of both the interaction and the Hamiltonians. 
Therefore, it is reasonable to say our model is quite universal in comparison with the previous 
models for quantum measurement. 

The system S plus the detector D forms a composite system C , the "universe" . In this sense, 
C is closed and thus its evolution can be described by an unitary operator U(t) governed by the 
total Hamiltonian 

H = H s + Hd + Hi=Ho + Hi 

By changing into the interaction picture and then backing to the original Schrodinger picture, a 
direct calculation gives the evolution operator 

u(t) = U0W1V); 

U I (t) = J2u n (t)\n)(n\ (2.4) 

n 

where 

M N 

U (t) = e^ Hat = 2_^e— —\n)(n\ ® [[ e s 

n=l fe=l 
N 

U n (t) = l[uW(t) (2.5) 

fe=i 

Vexp[±J*V n [x k (t)]dt}, for 0<t < r 
U [ *\t), for t>r 

Notice that V is the time-order operator and 

-iH k t iH k t 

Xk(t) = e i x fc e b (2.6) 

is the representative of the variables of -D in the interaction picture. 

It should be pointed out that the above mentioned evolution operator Ui(t) in the interaction 
picture just possesses the factorizable structure found in ref.[14, 15] by this author, which may result 
in von Neumann's wave packet collapse in the quantum measurement. 



UW(t) 



3. Quantum Decoherence at Finite Temperature 

If the system S were closed, the Schrodinger time evolution should not lead to the phenomenon 
of decoherence - wavefunction collapese defined by eq.(l.l). In fact, a pure state as a coherent 
superposition of some eigenstates of an observable of S can not evolve into a mixed state. This is 
because the unitary time evolution operator U(t) preserves the rank of density matrix for S, which, 
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however, has rank 1 for a pure state and has rank larger than 1 for a mixed state. We can also 
describe this impossibility in terms of the definition of quantum entropy in section ?. 

However, in our present model, S is considered as an open system interacting with D. Though 
the unitary evolution of the universe C formed by S plus D can not change C to a mixed state 
from a pure state, its induced effective evolution of S given by removing the variables of D may be 
non-unitary and thus the rank of the reduced density matrix can be changed in the evolution. In this 
sense, it is possible to realize the decoherence for quantum measurement in a quantum dynamical 
process. 

Assume that D is in a mixed state, e.g., an thermal equilibrium state at temperature T, when 
the interaction between S and D switches on. This initial state of D is denoted by the canonical 
distribution 

N 

Pd(0) = n®Y,Pn k \n k ><n k \ (3.1) 

fe=l n k 

where P nk is the classical probability of k'th particle of the detector in the states \n k > 

Let the initial state of S is a pure state, which is a coherent superposition 

|0>=^C„|n> (3.2) 

71=1 

of the energy eigenstates of S without the interaction with D. Then, we write down the density 
operator of the initial state for the universe C = S + D: 

p(0) = |$M®p D (0) = y] C m C*\m)(n\(g>p D (0) (3.3) 

Solving the von Neumann equation 

ih^p{t) = [p{t),H] (3.4) 
we obtain the density operator of C at time t 

P {t) = U(t)*p(0)U(t) 

Taking the partial trace of p(t) over the variables of D, we have the reduced density matrix 

M M 

Ps (t) = Tr D (p(t)) = VJ \C n \ 2 \n)(n\ + VJ (C n C*,F^ n (N,t)\n)(n \ + h-c) (3.5) 

n — 1 n>n 

where the off-diagonal terms are accompanied by factors 

-i(B n -E ,)t 

F n ,n'(N,t) = e r^Tr(U n ,\t)U n (t) PD (0)) 



iV 

nE < * it ^ i c)K >p « 



k=l rife 
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-»(E n -E ,)t " 

fe=i 

In the above expression, Tr^ means taking partial trace over the k'th variable of D. 

Now, we are trying to find a central condition, under which as the particle number N of D 
approaches infinity, F min (N,T,t) approaches zero for m ^ n, 

\F m .n(N,t)\^0 as N ^ oo (3.7) 

Eq.(3.7) determines the vanishing of the off-diagonal elements of p s (t) , that is to say, it will 
approximate the classical behavior of the open system 5* in the macroscopic limit with very large 
N. 

To this end, we use the eigenstates |a) of Hk(xk) with the corresponding eigenvalues s a (k) to 
rewrite the accompanying factor F n n > (N, t) = F . Then, 

n k n k 

Here, it has been taken into account that 

\(n k \U^{t)U^(t)\n k )\ 2 = l- Yl \{m k \U^(t)U^(t)\m k )\ 2 <l, fort = r (3.9) 

m k =£n k 

for an unitary operator U^)\t)Un\t) and n ^=n. In terms of the non-zero positive real number 

A n k n '(t) = -lhn\F^,(t)\ (3.10) 
the norm of accompanying factor F n n > (N, t) is expressed as 

Due to eq.(3.8), A™'" < 0; In general, under a reasonable condition that there is at least one 
non- vanishing off-diagonal clement in a given row or column of each factor U [ *}\t)Tj£ ] {t), A™'"(i)'s 
do not approach zero as k — > oo at t = r according to eq.(3.9). In this sense, X^fcLi A™'"(t) is a 
diverging series with the limit of infinity. In next section, some examples are present to obey the 
above mentioned condition explicitly. 

The above discussion shows the possibility of realizing the decoherence in a quantum dynamical 
process at the macroscopic limit with very large N. The above discussion does not depend on the 
specific forms of both the single particle Hamiltonian Hk(xk) and the interaction V n (xk)- It can 
be invoked to find a vast class of quite general dynamical models for quantum measurement to 
describe quantum decoherence - wave function collapse as the result of the dynamical evolution at 
macroscopic limit . 
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4. An concrete models with D made up of oscillators 



In this section, we will use a typical model to explicitly illustrate how the above approach works 
effectively for the dynamically- vanishing of the off-diagonal elements of the evolving density matrix, 
concretely speaking, under what kind of circumstance the series Y^kLi A™ n {t) m the exponential 
accompanying factor diverges into infinity. In this typical model, D is made up of N oscillators 
with Hamiltonian 

N 

H = ^fiuj k a+a k (4.1) 

i=l 

and the interaction between D and S is 

M N 
n—1 i— 1 

where the requirement fi m ^ fi n (for m ^ n) means that the coupling of D with S has different 
strengths for the different states \n> of S. 

In this sense, according to Wei-Norman's algebraic method [34], a factor Un\t) of the evolution 
matrix can be assumed as 

U£\t) = e /?(*) e A *(*)«ibe B *(*) * (4.3) 
Here, the coefficients F k (t),A k (t), and B k (t) to be determined satisfy a system of equations 

k k w ihuj 
ft{t)=BF{t)AUt) (4.4) 



It leads to 



and the real part of f£ (t) 



Al{t) = -Br{t) = ^[l-e^] (4.5) 



Mm)) = -\\m)\ 2 = -fja - cosc fe i) (4.6) 

is negative for uj k ^ nir/t. 

To master the kernel of the problem, we consider an simple case with zero temperature. The 
corresponding density matrix is pr>(0) = |0 >< 0|. Here, the ground state of D 

|0 >= |0i > <8|0 2 > (8i... (8 |0jv > 

is a direct product of the vacuum states |0j > (i = 1,2,---, N) of N oscillators. In this sense, 

= I < 0\UW(t)UW(t)\0 > | = e - A ™"W (4.7) 
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where 

2 -2 wjjt 

Ar(i) = 2(^ m - M») 2 (^ fc)2 2 = m(t)) (4.8) 
In the following discussion, we will detail the discussion about diverging of the scries J^k 
for various spectrum distributions of D. The most simple case is that D has a constant discrete 
spectrum, i.e., uok = u =constant,<jr/j = g =constant. In this case 



\F mn {N, t)\ = exp{-2N(p n - PnYj^ sin 2 — } (4.9) 

(nuj) z 2 

approaches zero as N — > oo except for the period points Lut = 2kir{k = 0, 1,2...) . Generally, the 
series J2k can repressed in terms of a un-specific spectrum distribution p{u)k) as 

^ = 2(, m -,„) 2 E PK ^ 2¥ (4-10) 

In fact, in the case of discrete spectrum, the spectrum distribution means a degeneracy that there 
are p(iVk) oscillators possessing the same frequency u>k- So long as the term 

pK)SfcSm 2 y/(ftwt) 2 

does not approach zero as k — > oo, the above series SR must diverge into infinity. For example, 
if LOk = kuj, p{iOk) oc k v+2 /gl,ri > 0,then this series must diverge into infinity for each term 
oc W sin 2 > 

Notice that f{M 1 t) = F m ^ n (N,0,t)l and means the complete coherence and the complete 
decoherence respectively. Because of the cut-off of frequency, the r\k {t) do not approach zero except 
for the period points t ^ ^"^^ as k — > oo. Also due to rjk > 0, the scries X)fcLi 7 7fe(*) mus t 
diverge into infinity, that is to say, f(N,t) — > according to eq.(3.4). Therefore, when the detector 
is macroscopic (N — > oo), the off-diagonal elements in ps(t) vanishes and the wavefunction collapse 
appears as the result of the dynamical evolution. It is realized from eq.(4.8) that, to realized the 
wavefunction collapse, we must constrict the interaction for each time measurement only to take 
place in a time interval r less than the oscillation period r c = Otherwise, the coherence terms 
suppressed by the factor fk(t) = F^ n (0, t) Owill be resumed at the common period points t = t c , 
Vk{tc) = and then \fk(t)\ = 1. The figure 1 with a constant spectrum ujk = u) shows how the 
decoherence appears in accompany with : \f(N,t)\ — > as N increases for t ^ t c , and how the 
coherence is resumed at t = i c . Figure 2 displays the same problems for the random spectrum 
distribution that the frequencies Uk take random value with cut-off. 

According the above analysis , one should let the interaction switch off before the coherence 
restores in the wavefunction collapse quantum dynamically. For this reason, the random spectral 
distribution and the constant spectral distribution are much appreciated for our present study for 
the QZE. For the case with continuous spectrum, some interesting circumstances can result from 
the concrete spectrum distributions. In the first example with 

p(u> k ) = 

9k 
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the series is convergence to a positive number proportional to time t 

sr = 5> r - 2(M % M " )2 f = ^-^ a * (4.ii) 



h z 7 ^ 2?i z 

This shows that the norm of the accompanying factor F mn (0,t) is an exponential decaying factor, 
i.e, 

\F m , n (0,t)\ ^e-"'""' 2 ' (4.12) 

This quite interesting result is much similar to that in the quantum dissipation [4] . As t — ► oo the 
off-diagonal elements of density matrix vanish simultaneously! 

In another example for continuous spectrum, the spectrum distribution is Ohmic type (e.g,see 
rcf.[4] ),i.e, 

N 2 

JH = S^%-Wj)=W (4-13) 



its alternative reformulation is 

*9k 9t 



(4.14) 



Then 

E\mn 2(/i m — /in) 2 S f°° . 2 Wfei K , 
A ™ — " T2 / sm ^- rfw fe = °°> f° rm ^ n (4.15) 

In summary, so long as we choose a suitable spectrum distribution of oscillators in the detector, 
the series ^ fc A™ n can diverge into infinity, that is to say, the dynamical evolution of S plus D can 
result in the wavefunction collapse or quantum decoherence, though the discussion in this section 
was proceeded with the oscillator detector, any detector (or environment) weakly coupling to system 
may be equivalent to a system of oscillators according to the proof given by Caldeira and Leggett 
[4]. Therefore, the discussion in this section do not lose the generality of the problem. 



5. Entropy Increment in Decoherence Process 

Since quantum decoherence decreases the information available to the observation about the 
quantum open system S, the quantum entropy 

S[p] = ~Tr{p\np) (5.1) 

as a functional of the density matrix p can be used to characterize the degree of decreasing infor- 
mation quantitatively. 

In comparison with the statistical thermodynamics, the decoherence process can be understood 
as an irreversible process in terms of the concept of entropy. In fact, if |A) is the eigenstatc of p 
with eigenvalue A, then eq.(5.1) is re-expressed as 

S\p] = ~ £ AlnA (5.2) 
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Obviously, the entropy is invariant under an unitary transformation and the time evolution of a 
closed system must not change its entropy. In our problem, the entropies of the initial and final 
states are zero and 

K M 

S[ P f]=- T J2\ C n\ 2 ^\Cn\ 2 (5.3) 
Z n=l 

respectively, manifesting that the decoherence process is certainly a process of increasing entropy. 
Notice that eq.(5.2) defines the maximum entropy of the system S for a given initial state (3.2) of 
S, which corresponds to complete decoherence. 

In a well-established theory for quantum decoherence, it is expected that the partially- and 
completely- decohering states are also described very well. For this reason, we calculate the cor- 
responding entropy of the intermediate state characterized by the reduced density matrix p s (t) for 
finite N. 

For the seek of simplicity, we first consider the two-state system with M — 2 physically. This 
case may be understood as spin — | precession interacting with a reservoir. In this sense, the reduced 
density matrix with finite N is explicitly written down 



p s = 



\c\\ 2 aczF 

C 2 C{F* |C 2 | 2 



(5.4) 



where 

F = F 12 {N 1 T) = e -£ ( =i A *W 
In terms of its eigenvalues Ai = |(1 + x) and A2 = 3(1 — a;), 



x = x(N) = y/l -4(1 - |F| 2 )|Ci| 2 |C 2 | 2 

= \]l - 4(1 - e~ 2 £" \ A *W)|Ci| 2 |C 2 | 2 
the explicit expression of entropy is 



(5.5) 



S(N) = + x)l n (l +x) + {\- x)l n {\ -x)- 2l n 2] (5.6) 



where 

IICil - IC2II <x < 1 



Since 



dS(N) K in (l+x) cQ 



dx 2(1- x) 

and x decreases as N increases, S = S(N) is a monotonically-increasing function in the above 
domain of x. When N = 00, x takes its maximum value so that there appears maximum entropy 

S ro = -y(|C 1 | 2 ln|C 1 | 2 + |C 2 | 2 ln|C 2 | 2 ) (5.7) 
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The above quantitative analysis shows us that the complete decoherence means the maximum of en- 
tropy. For the sufficiently large N, the average value A = TV -1 J2k=i Afc(i) i s roughly independents 
for N 

\F\=e- 2NA 

is approximated by exponential function of N. 

The explicit entropy functions for the partial decoherence in mesoscopic case with a finitely large 
N can not be easily solved for a general system of M(> 2) energy levels. However, for few cases 
with very large N or very small N, some approximation methods can be developed to calculate the 
entropy analytically. 

Foe the case with small N we define 

A mn (JV) = F mn (N) - 1. (5.8) 

Since the function F n ^ n = 1 and F m ^ n {N) approaches the unity for small N, the norm of A mi „(7V) 
is quite small in this case. Thus, the reduced density (3.5) can be decomposed into the unperturbed 
part 

Po = p(0) = Y^C m C*\m){n\ = \^\ 

and the perturbed part 

M 

Ap = &m,n(N)C m C* n \m)(n\ (5.9) 

m,n— 1 

The unperturbed part po denotes a pure state and then has M non-degenerate eigenstates \if>). 
Choose |^o) = to be the first one of basis vectors for the Hilbert space. Then, the other N- 
1 basis vectors \tpk) = P\k) (k = 1, 2, . . . , N — 1) to \ip ) may be constructed in terms of the 
complementary projection operator 

N 

P=l-|*)<*|= J2( 6 mn-C* m C n )\n)(m\ (5.10) 

m— 1 

It is not difficult to prove that such N — 1 vectors are linearly- independent. It follows from the 
Smite rule that N-l vectors 

\u k )=J2 S k'k\M, fc = l,2,..,AT-l 

k' 

are built as the orthogonomal basis for the complete space corresponding to the zero eigenvalues of 
Po- Therefore, the time independent perturbation theory determines the approximate eigenvalues 
of p = po + Ap up to of second order 

JV-l 

A ~l - l(*l A PK)| 2 + ... = 1-6+... (5.11) 

fe=i 

A fe = |<* \Ap\U k )\ 2 + ... = Sk + ■■■ ,k = 1,2, ... M — 1 
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Obviously, the normalization of density matrix still holds as X^feLo 1 ^ k = ^ un der this perturbation. 
A direct calculation from eq.(5.11) results in the entropy function in small N limit as 

S(N) = - 5) ln(l - S) + £ 4 ln^ 4] (5.12) 

fe k 

K M_1 

= - T e k*i a pk>i 2 
fe=i 

Another case that can be handled analytically by using approximation method is that with large 
N. In this sense, \F mn \ is so small that the off- diagonal parts of p (eq.(3.5)) can be regarded as a 
perturbation. The unperturbed part 

M 

p M = S ^\C m \ 2 \m)(m\ (5.13) 

m— 1 

represents the mixed state with the maximum entropy and |fc) and | C fe | 2 are just its eigenstates and 
the corresponding eigenvalues. Invoking the time-independent perturbation theory and regarding 

Ap = p-p M (5.14) 

as the perturbation, we write down the second's corrections for the eigenvalue 

AA [2] _ \ " \F mn \ 2 \C m \ 2 \C n \ 2 

AA " - \ \C m \ 2 ~\C n \ 2 (5 ' 15) 

for the first order solution A^ ' = |CV| 2 . Then, the approximate entropy is obtained up to of second 



order 

^) = -yE[l^l 2 ( ln l^l 2 + 2AA i 2) ( 5 - 16 ) 
1 k 



6. Application to Quantum Zeno effect 

May years ago the quantum Zeno effect (QZE) was theoretically proposed by Misra and Sudar- 
shan [25] based on the postulate of wavefunction collapse. It is argued that an unstable particle 
will never be found to decay when it is continuously observed, more generally speaking, a frequent 
measurement inhibits the transitions between quantum states. Recently, Itano, Heinzen, Bollinger 
and Wineland (IHBW) [26] , have reported that they have observed the QZE in an experiment about 
atomic transition based on Cook's proposal. They claimed that the freezing of stimulated transition 
probability appeared when the two-level atomic system is subjected to frequent measurements of the 
population of a level. Then, studies of the QZE have attracted much interest over last years[27-37]. 
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Among these discussions, Petrosky,Tusaki and Prigogin's (PTP's) remarkable work[27-28] showed 
that the result in IHBW'S experiment can be recovered through conventional quantum mechanics 
and do not involve a repeated collapse of wavefunction. This is quite similar to Peres's observa- 
tion[37] that the a modified Hamiltonian may mimic the wavefunction collapse to slow down the 
quantum transition and then realized the QZE in a pure framework of quantum mechanics. It 
seems that the conclusion of IHBW's experiment is challenged by the theoretical analysis of PTP's 
work, and the same result can be obtained in term of both the use of "wave- function collapse" by 
frequent measurement , or the use of the Schrodinger's evolution as a pure quantum dynamical 
process [29,30]. Thus, It is difficult to say whether the IHBW's experiment has proved the existence 
of the QZE or not. 

To compromise the above mentioned different standpoints in a reasonable framework, we will 
reconsider the QZE and its corresponding IHBW's experiment from a distinct point of view. Since 
the wavefunction collapse for quantum measurement can be regarded as a quantum dynamical pro- 
cess in certain circumstance according to the above discussion, it is natural to understand the QZE 
and IHBW's experiment as the results of a measurement monitoring the system continuously, but 
the wavefunction collapse characterizing measurement can be realized here quantum dynamically. 

6.1 Model Hamiltonian for QZE and its Evolution Matrix 

Now, we present a dynamical model for the QZE. Let the measured object be a two-level systems 
S with Hamiltonian 

if = §(|i><2| + |2xrj) (6.1) 

where |1 > and |2 > are the ground and excited states respectively. SI is Rabi frequency of the 
external field coupling with atom. The continuous measurement for a given time interval T is imaged 
as a limit of the L times successive measurements at times t = kT/L (0 < k < L) with TV — ► oo . 
For fc'th measurement, the interaction of measuring instrument-detector D on the system turns on 
at time tk = kT/L and then turns off at at time tk + t for k = 0, 1, 2, ... . Here, the r is very short 
in comparison with T/L. This case is quite similar to the remarkable experiment by Itano et.al [26]. 
In their experiment involving three levels of 9 Be + ions, a on-resonance radiation frequency field is 
applied to 9 Be + with L shot on-resonance optical pulses-measurement pulses to perform a quantum 
measurement. Each optical pulse results in the off-diagonal elements of the density matrix to zero 
with the postulate of collapse wavefunction. 

In our model, the instrument D consists of N harmonic oscillators with free Hamiltonian 

JV 

H D = ^ huiajai (6.2) 
where and a.j are creation and annihilation operators for z's mode of boson states respectively. 
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To realized wavefunction collapse, the frequency cut-off must be introduced for the spectral distri- 
bution of the harmonic oscillators, that is to say, there exists a up-boundary oj c for the frequencies 
{u)\,u)2, un}. The model-interaction 

1 N n 

Hi = g + >< !| - I 2 >< 2 I) cos ^ * ^ ~ sin >< 1| - I 1 >< 2 I) + 1] 

is turned-on and turned-off by the switching function 




9, if t k <t < t k + t 
0, otherwhere 



Based on the same method as that in section 3, we obtain the exact solution U(t) for the 
Schrodinger equation of the total system C in terms of the free evolution matrix 



N 



e 



k=i 



cos^t -isin^t 
-isin^t cos^t 



(6.4) 



U (t) = e -*E< a ^ 4 • e'^ = U D ■ U s = (J[ , 

and the interacting evolution matrix 

U e (t) = U e (t)\l >< 1| + \2 >< 2\U e (t) = J[ e F ^e A ^ a+ e Bk ^ a = \[ U [k] (t) (6.5) 

where 



A I 



k=l 



M 



fe=l 



A k {t) = -B* k (t) = J* i^e Wkt dt = J^E e ™ ktdt + 

iwTl 

= ^-[(1 ~ e^ kt )e^ + (1 - e^ kT ) ~ & 
nio 1 — e i 

F fc (t)=BJE(t)A fc (t) 



(ft] 



(6.6) 
(6.7) 



for U < t < U+T. Here, we have used the properties of the switching function g(t). If r << 1, then 
the approximate function Ak(t) is 



A fc (i) = |-(l-e^V 



It follows from the above equations that 



Vk(t))Re 



[Fk(t)] = - l -j\X k {t)Ai{t)+A k {tyA k {t 



))dt = -\\A k (t)\ 



(1) 



(6.9) 



Notice that real part of F k (t) is negative and e Re ^ Fk ^ is not larger than unity. 
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6.2 QZE based on Dynamical Collapse 



For the measured system S with an initial pure state 

\ip >= d\l > +c 2 |2 > 

and the measuring instrument D with initial mixture state described by a density matrix 

N 

pd=yi po{k), 
fe=i 

the initial state of the composite system C formed by S plus D is a mixture state with density 
matrix 

p(0) = |V> >< VI ®Pd 

Then, the state of C at t(< t\) is 

p(t) = U (t)U e (t)p {0) U+(t)U+(t) 

Now, we understand the process of measurement as a procedure to determine the reduced density 
matrix, which contains the total information of the measured system S. The first measurement 
results in a reduced density matrix for S 

Ps = Tr D p(t) = U s (t)(\ Cl \ 2 \l X 1| + |c 2 | 2 |2 >< 2\ + f(N,t)\l X 2| + f(N,t)*\2 >< l\)U+(t) 

(6.10) 

by taking the trace for the variables of the detector D where the coherence factor 

N N 

f{N,t)=Tr D {U e {t) PD ) = Y[f k (N,t) =l[Tr k [U^(t)p D (k)} 

k=l k=l 

is factorizable. This factorization is crucial for the appearance of wavefunction collapse . In terms of 
the probabilities P n (k) of its distribution on the Fock state \n > for the k'th oscillator in measuring 
instrument, po(k) can be re-expressed as 

Pn{k) = ^P n {k)\n >< n\ 

we prove 

\f k (t)\ - \Tr k {U^{t) PD {k))\ = \Y,Pn{k) < n\UW(t)\n > | < £P„(fc)l < n\U (k \t)\n > | < = 1 

with a similar discussion to that in section 3. Thus, the N- multiple product of all |/fc(i)|'s must 
approach zero as N — > oo unless most of |/fc(i)|'s are unity simultaneously. It will be illustrated 
by the following typical example that the case that all or most of |/fc(t)|'s are unity is rather special 
and this case can be eliminated by choosing a spectrum distribution with cut-off of frequency. In 
this example, we take po(k) to a pure state |0 >< 0|, obtaining 

f(N,t) =< 0\U T (t)\0 >= e - ^! (6.12) 
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where |0 >= |0i > ® | O2 > ® ■•■ ® |0m > is the vacuum state of the detector D consisting of M 
oscillators. 

1 o 2 

n k (t) = -ite[F fc (t)] = -|A fe | 2 = (i _ C0SLJfc t) ( 6 .i3) 

Notice that eq.(6.13) is quite similar to eq.(4.8). Then,, when the detector is macroscopic 
(N — ^ 00), the off-diagonal elements in ps(t) vanishes and the wavefunction collapse appears as the 
result of the dynamical evolution, that is 

P(t) ^ps = U s (t)(\c\ 2 \l >< 1| + \c 2 \ 2 \2 >< 2\)U+(t) (6.14) 

In the following discussion, we use the above developed theory for the quantum dynamical model 
of wavefunction collapse to recover the QZE. 

According to eq.(3.4), after the first time of measurement, the density matrix is 

Ps(h = j) = >< 1| +4 21 |2 >< 2| + AW|1 >< 2\ + A^*\2 >< 1|) + oW(l) (6.15) 



where 



C W = | Cl | V + | C2 | 2 /3 2 C W = | C2 | V + | Cl | 2 /3 2 AW = ^| Cl | 2 /3 - l -\c 2 \ 2 (6.16) 

flT flT 
a = a(L) = cos—, f3 = (3(L) = sin— 

Notice that the last term in eq.(4.1) is 

[ K ] (jj) = U^[f(N,t)\l >< 2\ + f*(N,t)\2 >< l\]U+(t) (6.17) 

which vanishes as N — ► 00, the first term in eq.(3.1) 

C/ s (* 1 )(|c 1 | 2 |lxl| + |c 2 | 2 |lxl|)C/+(t 1 ) 

represents the mixed state with complete decoherence. Subsequently, the measurement of the first 
time cancels the off-diagonal term of p s (Q) — X VI as N — ► 00. The measurement of second 
time will cancel the off-diagonal terms A^\l X 2| and ^4M*|2 X 1| in eq(3.1). Similarly, after 
the measurement of fc'th time 

Ps [k] = |4 fc] | 2 |l >< 1| + |c 2 fc] | 2 |2 x 2| + A^\l >< 2\ + A^*\2 x 1| +0[ 1 ](1) 

= t/5(t fc )(| C i ''I'll X 1| + \c [ t 1] | 2 |2 x 2\)U s {t k ) + + OW(-L) (6.18) 
where c^' and c 2 fe ' arc determined by the recurrent relations 

[k] [fe-1] 2 , [fe-l] 2 [k] [fe-1] 2 , \k-l] a 2 , c1 m 

q J = q a + c 2 p,c 2 = c 2 a + c{ p (6.19) 

where 

= U,{t k -x)\A^f{N,t k )\l X 2\+A^*f(N,t k )*\2 >< l|]tf+(t fc -i) 
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also disappear as N — ► oo 

It is quite difficult to solve J-^ and explicitly from cq. (6.19), but wc can invoke the computer 
simulation to evaluate the variations of and as the measurement times L increase. If the 
system is initially in pure state |1 >, the initial conditions are C\ = 1 and c 2 = 0. It is illustrated in 
Figure 3 that , as L — ► oo, the distributions |c^ L '| — ► oo while \c^\ — > 0. Notice that |c^ L '| first 
decrease and then increase to 1. This means that, for a macroscopic instrument (N — > oo), the 
system will be forced back to the state 1 1 > as the successive measurements become continuous for 
L — > oo. This just the QZE! 



6.3 "Transition" of Quantum Entropy for QZE 

Now, we consider the entropy of the measured system S as a function of the times N of mea- 
surement. Considering that the entropy is invariant under an unitary transformation due to the 
original definition , one can obtain the entropy for the QZE after L'th measurement 

S(L) = S[p s (L)\ = -^(\c[ L] \\n\c[ L] \ + |4 LI |In|4 L1 l) (6-20) 

The entropy S = S[ps(L + 1)] takes its maximum value 

S max {L + l) = ^ln2 (6.21) 

only when 

i4 l] m4 l] i. 

According to the recurrent relations eq.(6.19), the above equation leads to 

\cf\ \£\ = (a 2 - \cf-% = (a 2 - f )\c\ cl) 

Then, it is easily observed that the above equation holds only when 

a 2 (L)=(3(L) 2 

or 

^ = |(2i + l),/ = 0,±l,±2,... (6.22) 
The entropy S takes several maximum values 

S m (L) - ~[cf ] lnc[ L] + cf \ncf] (6.23) 

where 

£=[^+1^]. l = 0,±l,±2,.... (6.24) 
Here, [x] denotes the integer part of the real number for x > 0, if x < then [x] = 0. 



20 



Let us consider some special examples for the above discussions. If uT = n, L = 2 ( for I = 0), 
and there only one point for the maximum entropy 

5 L+1 = -^(|4 2 ]| 2 ln| C [ 2l | 2 + |4 il | 2 ln|4 il | 2 ) (6.25) 

Generally, if wT = kn for positives integer k , 

^ r 2/c n ,.2/c. r 2/c. r 2/c n .„ 

L = lWTi ] = 2k ^^"^2xTi ] (6 ' 26) 

where X < k — 1, that is to say, there are fc points 2k, [^], [ 2 2 +i ] f° r the maximum values 

of entropy. 

For the general case with uoT = kT, the above analysis shows us that, if L is less than the 
critical value L c = 2k, the variation of S(L) is a "random function" of L, which is not monotonic; 
however, when L is larger than the critical value 2k, S(L) is a monotonically-decreasing function of 
L . When ujT is not an integer times of n ,the critical point for L is ^1+1) 7r ' ^ ucn a feature of 
"transition" from random to regularity for quantum entropy in the QZE is illustrated in Figures 4 
and 5. Physically, the QZE defines a transition of the information entropy from random to regularity. 



7 Generalization of Cini's Model for Quantum Measurement 

It has to be pointed out that the correlation between the states of the measured system and 
that of the detector has not been emphasized well in the original Hepp-Colcman model and its 
generalization . This problem was well analysed by Cini with a beautiful dynamical model [39]. 
The present investigation is to emphasize on both the wavefunction collapse and the correlation 
collapse for a generalization of the Cini model. In fact, the correlation between the states of 
measured system and the detector is crucial for a realistic process of measurement, which enjoys a 
scheme using the macroscopic counting number of the measuring instrument-detector to manifest 
the microscopic state of the measured system. 

The original Cini model for the correlation between the states of measured system S and the 
measuring instrument-detector D is build for a two- level system interacting with the detector D, 
which consists of indistinguishable particle with two possible states luq and u>i . For the two states 
u + and U-, the detector has different strengths of interaction with them. Then, the large number 
N of "ionized" particle in the ionized state u>i transiting from the un- ionized state luq shows this 
correlations. 

7.1 Generalized Cini Model 



21 



In this section, we wish to generalize the Cini's model for the M-level system. The measured 
system S with M-levels has the model Hamiltonian 

M 

H s =Y / E k\®k><<S>k\ (7.1) 
k=l 

Where |<&fc > are the eigenstates corresponding to the eigenvalues Ek (k = 1,2, ■■■,M). The 
detector D is a two-boson-state system with the free Hamiltonian 

Hd = huiia^ai + hw 2 a 2 a 2 (7.2) 

where ai and a/" are the creation and annihilation operators and they satisfy 

[cLi,af] = S i:j , [a u aj] = [a+, a+] = (7.3) 

In the Schrddinger representation, the interaction is described by 

Hi{t) = Y J 9e~ r,t {wn\$n >< ^ n \)(e l{ ^- UJl)t a+a 2 + e^-^a+a^ (7.4) 

n 

where the non-degenerate weights w n represent the different strengths for the different states |$„ > 
of the system. The exponential decay factor e~ vt for r\ > is here introduced to turn off the 
interaction after suitable time so that the coherence can not restore in the evolution process. This 
point can be explicitly seen in the following discussion. The introduction of time-dependent factors 
e ±i(oJi-w 2 )t j s q m te similar to that in ref.[13] where these factors are used to describe the energy 
exchange due to the presence of the free Hamiltonian Hp- Notice that there was not the free 
Hamiltonian for detector in the original Hcpp-Colcman model such as Hd in our present model. 
Let H = Hs+Hp+Hi be the total Hamiltonian in the Schrddinger representation for the composite 
system formed by S plus D. Transforming the problem into the interaction representation with the 
evolution operator 

U (t) = exp[^-(H s + H D )t] 
in 

one has the interaction potential 

M 

V I (t) = e- r > t gJ2W n \$ n ><$„|(a+a 2 + a lfl +) (7.5) 

ra=l 

In order to diagonalize Vi(t), we invoke the canonical transformation as in ref.[5] 

ai = -J=(6i -b 2 ), a 2 = -J= {h + b 2 ) (7.6) 

where the new boson operators bi and 6+ satisfy the same bosonic commutation relation. In terms 
of these operators, Vj (t) is rewritten as the diagonal form 

M 

V T (t) = ge-^ Wn\$>n >< $ n \(bth - b+b 2 ) (7.7) 

n=l 
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Then, considering that the interaction part Vi(t) commutes with each other at different time, ie. 

[V I (t),V I (t')] = 

one can express the evolution operator 

If* M 
Ui(t) = exp[- / Vt^dt'] = V expi-it^gWnib+h - b+b 2 )]\$ n >< $ n | 

m J o „=i 

M 

= ^(7„(t)|$„ ><$„| (7.8) 



where 

1 - e-rt 



n 



takes the real time as r\ — > 0. It can be regarded as the 77- deformation of time t; when t — > 00, t,, — > ^. 



7.2 Correlation of states from Evolution of State 

Now, we consider the evolution of the total system starting with an initial state at t=0 

M 

|*(0) >=^C fe |$ fe >®\N,0> (7.9) 
fc=i 

where 

\m,n>=^ T \0> (7.10) 
v mini 

denotes a Fock state of two-boson system which denotes that there are n particles in " ionized" state 
and m particles in un-ionized state. It is hoped to manifest a correlation between the states \m, n > 
of detector and the state |$; > of the system in a dynamical evolution of the state \^(t) > in some 
limiting case so that one can read out the state 1$; > from the manifestation of the state |m, n > 
of the detector. Notice that the eigenstates of the operator 



6 = b+h - b+b 2 (7.11) 



are 



I A , N - A} = , 1 6+ A 6+ (JV " A) lO > 
y/M(N- A)! 1 2 

= — > > —= V V ' =m,n> 7.12 
V2 ^ V(A - m)\(m\)(N - A - n)!(n!) 

with the eigenvalues 

£Ar(A) = 2A-Af (7.13) 

where A = 0, 1, 2, • • • , TV. for a given integer N. The original Fock state can be expended in terms 
of |A,AT-A> as 

1 N N \/W(-i\ N - x 

|JV,0 >=(-=) Y ,' { ' =\X,N-X} (7.14) 
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Then, one obtains the wavefunction at time t 



|*i(t) >= C/}(t)|*(0) > 

JV M 
A=0fe=l VZ 

/ M(-i) w - A 



e -iflt,^(2A-JV)| $fe > ^N- A > (7.15) 



Then, one has 



where 



\/A!(iV - A)! 

M 

|*j(t) -^C fe |$ fc >®^a„(t,fc)|n,iV-n> (7.16) 



M N 



k=l n=0 



a n (t, k) = \ V = cos n {gW k t ri )sin N - n {gW k t v ) (7.17) 
\Jn\(N — n)\ 

Obviously, the probability ol finding n "ionized" particles in second bosonic mode is 

P n = Kl 2 = n , { n- n) ! ^ (dWkt^sin 2 ^-^ (gW k t„ ) (7.18) 

or 

P n = C%p n (l-p) N - n (7.19) 

where 

when N is very large so that the Stirling formula is valid, it can be proved that when n k — n k — 
Np k (t), the probability p n has its maximum. 

Pn k = Cg(^)(^±) N -^ (7.20) 

Notice that the derivation is the same as that in ref. [39] , but n k depends on the index k. As proved in 
rel.[39], P nk is very strongly peaked around its maximum Pn k , which becomes unity when N — > oo. 
Therefore, if the detector is very macroscopic (N — > oo), then P nk (N — > oo) = <5„ fcfifc that leads to 

M 

|*/(t) >w^oo= E Cfc l* fe > A?" - fife(t) > (7-21) 

fe=i 

When fih(t) ^ nk'(t) for fc ^ fc', a one-to-one correlation between the states |$fc > of S and the 
states of D. In this sense, if the detector is found in the state \fik(t), N — fik(i) >, it can be concluded 
that the system is in this state \<& k > A realistic detector must have a good fringe visibility, which 
can manifests the macroscopic differences between any two states of \fik, N — n k > for different k. 
It required that there is not the considerable overlap between \n k , N — h k > and \n k i , N — n k > > for 
k 7^ k'. In fact, if n k = fik' Then 

g W k t v = gWk't n + nw, n = 0, 1, 2, • • • 
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or 

gt v {W k -W k >) = nir 

Otherwise, if one lets the interaction between S and D decay very fast so that the limiting time 

I< 1 

v ~ g{w k -w k ,) 

for any k ^ k' , then fifc(i) ^ fik'(t) for any long time evolution. Notice that the above mentioned 
problem of overlap of correlation states is a disadvantage in the original Cini's model, where there 
is not the decay factor of interaction. So, in the time t n = the states \n,N — n > and \N,0 > 
are completely overlap. Here, it is the case of two levels with W2 = W = -j=,W\ = 0, thus, at 
t = t n the correlations vanish for the original Cini's model. Introducing the decay factor e~ vt in the 
interaction is the key point to avoid vanishing of correlation in our model. In fact, such decay of 
interaction can appear in realistic physics. For example, an atom is prepared in a microwave cavity 
loaded with an electromagnetic field which can decay at a suitable rate. In this example, the atom 
and the cavity are regarded as the system and the detector respectively. Notice that this example 
is quite useful for the studies of atomic cooling [40] . 

7.3 Wavefunction Collapse in Cini Model 

We can also use the above generalized Cini's model to describe the wavefunction collapse for 
the M-level system quantum mechanically. Let the system D be initially prepared in a coherent 
superposition of M-level 

M 

|*>=]Tc fe |$ fe > (7.22) 
k=i 

and the detector be adjusted in the initial state \N,0 >, then the density matrix for the initial state 
of total system is expressed as 

po = CkC *k' \® k >< ® k ' \®\ N >° >< N > °l ( 7 - 23 ) 

k,k> 

Using the evolution operator Uj(t) in the interaction representation, one formally write down the 
density matrix for the total system at t 

p(t) = U!{t) P{ Q)Ut{t) 

Because we are only interested in the final state of the system other than that of the detector for 
the consideration of wavefunction collapse, so we must take trace for the variable of detector in the 
total density matrix to obtain a reduced density matrix for S 

p s (t) =Tr D p{t) = <m,m'\p(t)\m,m / > 



25 



= Y^\Ck\ 2 \^ k ><<M 
k 

+ £ e- u ^ E "- E ^l h cos N [9{W k - W k ,)t v }C k C* k ,\$ k >< $ fc ,| (7.24) 
where we have used 

N 

^K(t,fc)| 2 = l (7.25) 

n=0 

and 

N 

J2an(t,k)a* n (t,k') = cos N [g(W k -W kl )t v ] (7.26) 

n=0 

Notice that each off-diagonal element in the density matrix is accompanied by a time-dependent 
factor 

F N (k,k') = cos N [g(W k - W k ,)t n ] (7.27) 

which is a N-multiple product of factors cos[g(W k — W k ')t v ]. Recalling that due to the existence of 
the strong decay factor e~ vt so that 

holds for any k ^ k' , we observed that the deformed time t n changes from t v = to t v = | as the 
real time changes from t = to t — ► oo respectively. In this sense, 

9(W k -W k ,)t v <ir (7.28) 

and then 

0< \cos[g(W k - W fe 0g| < 1 

or 

\cos[g(W k - Wy)^ = e-^\f k (t v ) > (7.29) 

This observation leads to 

F N {k,k') = e- N f^ 

which obviously approach zero as N — > oo. Therefore, when the detector is macroscopic, (N — > oo) 
the off-diagonal terms of the density matrix vanish and the wavcfunction collapse is realized quan- 
tum dynamically. 

To end this paper, we present some comments on the above discussions. Though this paper 
provides one with an extensive generalization and the unified description for a number of dynamical 
models ( e.g., the HC model ) for the quantum dccohcrcncc, we have to say that a disadvantage 
in the original model still exists in the present models. This is the oscillation of ul!^(t) U^f may 
enable most of the factors F}® ,(T,t) in eq.(3.6) to become unity at a specific time t = To and 
thus the whole accompanying factor can not approach zero at this time. To suppress such kind of 
oscillation so that the decoherence appears in dynamical evolution, a phenomcnological method is 
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to use the switching function g(t) in the interaction (2.3). However, the microscopic mechanism of 
this switching is not clear for us. We believe that the quantum dissipation caused by the detector 
or environment is a possible way to introduce such a switching mechanism microscopically. For the 
concrete example that D is made up of harmonic oscillators, the discussion in section 4 showed 
that this kind of dissipation may result from the specific distribution. How to realize the quantum 
decoherence directly through quantum dissipation is still an open question. 
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Figure Captions 



Figure 1: The norm \f(N, t)\ of coherence factor as the function of time t for the case with spectral 
distribution uik = u- Figures 1-a, 1-b, and 1-c correspond to the times of measurement equal to 5,10 
and 500 respectively. [Here , u k = u = 0.1, jj-^j = 0.00001]. For larger N (e.g. iV=500), the coher- 
ence enjoyed by \ f(N, t)\ almost disappear for ih <t < tfc+i- However, at th — ^j 1 , (k — 1, 2, N) 
the coherence is resumed in a very shot time. 

Figure 2. The normal of coherence factor for the spectral distribution , the uik is random with 
a cut-off frequence uj c =. The time for resuming the coherence is t c (k) = ^jjjz. 

Figure 3. The probabilities |c^| and \c^\ finding the system at states |1 > and |2 > respec- 
tively. 

Figure 4. The entropy as the function of times L of measurement for uoT = 3tt, it represent a 
"transition" of the information entropy from the random to regularity. Figure 4-b is part of Figure 

4- a. 

Figure 5. The entropy as the function of times L of measurement for ujT — 10n, it represent a 
"transition" of the information entropy from the random to regularity. Figure 5-b is part of Figure 

5- a. uoT = 107T, and Figure 5-c,d with uiT = lOn. 



28 



References 



1. J. von Neumann, Mathemstische Gruandlage de Quantumechanik, Berlin: Springer, 1993 

2. V.B.Braginski and F.Y.Khalili, Quantum Measurement, Cambridge University Presses, 1992 

3. Tso-hsiu Ho, Physics (in Chinese), 22(1992) 419 

4. W.H.Zurck, Phys.Today, 1991,No.lO, p 36. 

5. B. Kramer, Quantum Coherence in Mesoscopic Systems, New York: Plenum, 1991. 

6. M.Gell-Mann and J.B.Hartle, Phys.Rev.D, 47(1993), 3345 

7. A. J.Leggett, S. Chakravarty, A.T. Dorsey, M.P.A. Fisher, A. Gary and W.Zwerger, Rev.Mod.Phys., 
59(1987), 1. 

8. J.P.Paz, S.Habib and W.H.Zurck, Phys.Rev.D., 47(1993),488 

9. L.H.Yu and C.P.Sun, Phys.Rev.A, 49(1994), 592 

10. C.P.Sun and L.H.Yu, submitted to Phys.Rev.A. 

11. K.Hcpp, Hcv.Phys.Acta, 45(1972), 237. 

12. J.S.Bell, Hcv.Phys.Acta,48(1975), 93. 

13. H.Nakazato and S.Pascazo, Phys.Rcv.Lett., 70(1993), 1 

14. C.P.Sun, Phys.Rev.A, 48(1993),878 

15. C.P.Sun, Chin.J.Phys., 32(1994), 7 

16. N.G.van Kampcn, Physica A,153(1988), 97 

17. M.Namik,S.Pascazio, Phys.Rev.A44.(1991),39. 

18. M.Namik, Found. Phys.Lctt,5(1992) 265 

19. B.Gaveau and L.S.Schulman, J.Stat.Phys. 58(1990), 1209 

20. C.P.Sun,J.Phys.A,21(1988)1595. 

21. C.P.Sun, Phys.Rcv.D.38(1988),2908. 

22. C.P.Sun, High Energy Phys.Nucl.Phys., 12(1988) 252 

23. C.P.Sun,Phys.Rcv.D,41(1990),1318. 



29 



24. C.P.Sun, Physica Scripta,42(1993)393 

25. . B.Misra and E. Sudarshan,J.Math.Phys. 18(1977)756. 

26. . W.M.Itano, O.J.Heinzen, J.J.Bollinger and D.J.Wineland, Phys.Rev.A 41(1990)2295. 

27. . T.PctrovskyS.Tasaki and I.Prigoginc. Phys.Lett. A 151(1990)109. 

28. . T.Petrovsky, S.Tasaki and I.Prigoginc, Phyica A 170(1991) 306. 

29. . S.Inagaki, M.Namiki and T.Tajin, Phys. Lett. A 166(1992)5-12. 

30. . S.Pascazio, M.Namiki, G.Badarek, H.Ranch, Phys. Lett. A 179(1993)155. 

31. . D.Homc,M.Whitakc, Phys. Lett. A 173(1993)327. 

32. . R.Onofrio,C.Presilla and U.Tambini, Phys. Lett. A 183(1993)135. 

33. . T.P.Altenmiiller and A. Schenzle, Phys.Rev. A 48(1993)70. 

34. . M.J.Gagcn, H.Wiscmcn, G.J.Milbarn, Phys.Rev. A 48(1993)132. 

35. . T.P.Altenmiiller and A.Schcnzlc, Phys.Rev.A 49(1994)2016. 

36. . G.S.Agarwal and S.P.Tewari, Phys.Lett. A 185(1994)139. 

37. . A.Pcrcs, Am.J.Phys.48(1980)931. 

38. . J.Wci and E.Norman, J.Math.Phys., 4A(1963) 5 

39. . M.Cini, Nuovo Cimcnto,73B(1983),27. 

40. . T.Zaugg, M,Wikcns, P.Mcystre and G.Lcnz Optics Commun. 97(1993)189. 



30 



